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ABSTRACT 



We report the detection of emission in the 7=1-0 rotational transition of hydrogen fluoride (HF), together with observations of the 
7 = 1-0 to 7 = 3- 2 rotational lines of H 35 C1 and H 37 C1, towards the envelope of the carbon star IRC +10216. High-sensitivity, 
high-spectral resolution observations have been carried out with the HIFI instrument on board Herschel, allowing us to resolve the 
line profiles and providing insights into the spatial distribution of the emission. Our interpretation of the observations, with the use of 
radiative transfer calculations, indicates that both HF and HC1 are formed in the inner regions of the envelope close to the AGB star. 
Thermochemical equilibrium calculations predict HF and HC1 to be the major reservoirs of fluorine and chlorine in the atmospheres 
of AGB stars. The abundances relative to H 2 derived for HF and HC1, 8 x 10~ 9 and 10~ 7 respectively, are substantially lower than 
those predicted by thermochemical equilibrium, indicating that F and CI are likely affected by significant depletion onto dust grains, 
although some chlorine may be in the form of atomic CI. The H 35 C1/H 37 C1 abundance ratio is 3.3 ± 0.3. The low abundance derived 
for HF in IRC +10216 makes it likely that the fluorine abundance is not enhanced over the solar value by nucleosynthesis in the AGB 
star, although this conclusion may not be robust because the HF abundance we derive is a lower limit to the elemental abundance of 
F. These observations suggest that both HF and HC1 should be detectable through low 7 rotational transitions in other evolved stars. 
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1. Introduction 



identification — molecular processes — stars: AGB and post-AGB — circumstellar matter — 



Light hydrides such as HF and HC1 are very difficult to study in 
space by means of their rotational spectra, which lie in the sub- 
millimeter and far-infrared spectral regions that are hard if not 
impossible to observe from the ground owing to the opacity of 
the terrestrial atmosphere. Moreover, the prospects for detection 
are limited by the low cosmic abundance of fluorine and chlo- 
rine, although HF and HC1 are expected to be the major reser- 
voirs of each of these elements in many astrophysical environ- 
ments. 

Since it is one of the lightest molecules, HF has a rotational 
spectrum that can only be observed from space. It was first de- 
tected with the Infrared Space Observatory in Sagittarius B2 
dNeufeld et al.l 1997). Observations of HC1 have been conducted 
in molecular clouds dBlake et al.l 1985: IZmuidzinas et al.l 1995: 
ISchilke et ail 1995; ISalez et al.l 1996) through its fundamental 
rotational transition, lying at 625.9 GHz and still accessible from 
ground, and in the diffuse cloud £ Ophiuchi (Fede rman et al.l 
1995) through an electronic transition at ultraviolet wavelengths 
using the Hubble Space Telescope. That HF and HC1 are very 
stable molecules makes them abundant in the atmospheres of 
cool stars and in sunspots, where they have long been ob- 



served through ro-vibrational transitions in the infrared region 
flHall & Noyes| 1969; |Hall & Noyes] 1972; |Ridgway et aL] 1984; 
IJorissen et al.l 1992). 

With the launch of the Herschel Space Observatory, it has 
been possible to perform very sensitive observations of the 
pure rotational lines of HF and HC1 in different types of in- 
terstellar clouds (se e INeufeld et all 2010; |Phillips et al.| 2010; 
ISonnentrucker et aXI 2010; ICernicharo et al.l 2010a: ILis et al.l 
2010). In this Letter, we present sensitive observations of the 
7=1-0 rotational transition of HF and of the 7 = 1 - to 
7 = 3-2 transitions of H 35 C1 and H 37 C1 towards the carbon star 
envelope IRC +10216. This is the first time that HF and three 
lines of the two major isotopologues of HC1 are observed at high 
spectral resolution in IRC +1021(fl allowing us to resolve the 
line profiles and gain insight into the spatial distribution of the 
molecules. 



2. Observations 

Our HIFI observations of the 7=1-0 line of HF and of the 
7=1-0 to 7 = 3- 2 lines of H 35 C1 and H 37 C1 were obtained 
in May and from October to December 2010, in the context of 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



1 Several transitions of HC1 were observed at low spectral resolution 
with the SPIRE and PACS instruments JCernicharo et al.l 2010b ), and 
the 7=1-0 line was observed with the CSO telescope ( Peng et al. 2010). 
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Fig.l. HIFI spectrum of IRC +10216 showing the 7 = 1-0 
rotational line of HF (black histogram; spectral resolution is 1.5 
MHz) and the line profile resulting from the radiative transfer 
model (blue line). 

GT1 programmes to perform a line survey and to search for 
light hydrides in IRC +10216. Data were taken in double beam- 
switching mode with a spectral resolution of 1 . 1 MHz, and pro- 
cessed using the standard Herschel pipeline up to level 2, which 
provides fully calibrated spectra. The local oscillator was shifted 
in frequency to identify any emission from the image band. Total 
integration times were 43 minutes for the HF line, and ranged 
from 10 to 21 minutes for the HC1 lines. The final spectra were 
smoothed to a spectral resolution of 1.5 MHz and have antenna 
temperature rms noise levels in the range 0.04-0.1 K. For details 
about the data reduction, we refer to Cernicharo et al. I (20 10c). 



3. Results 

The 7 = 1 - line of HF and the 7 = 1 - to 7 = 3 - 2 lines of 
H 35 C1 and H 37 C1 are all clearly detected, free of contamination 
by other lines (see Fig. Q] and Fig. 01. The differences between 
the observed and calculated frequencies are not negligible, up 
to ~10 MHz (see TableQ}. They are, however, within the errors 
and, expressed in Doppler equivalent velocity, have values in the 
range 1.3-1.9 km s _1 , thus are much smaller than the observed 
line widths (21-31 km s" 1 ). The detection of HF, H 35 C1, and 
H 37 C1 in IRC +10216 is therefore quite secure. 

The HC1 hyperfine structure due to the nuclear quadrupole 
of the 35 C1 and 37 C1 nuclei can only be barely distinguished in 
the 7=1-0 lines, but not in the higher 7 lines owing to the 



Table 1. HF and HC1 observed line parameters in IRC +10216 
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and Cazzoli & Puzzarini (2004), respectively. b v cxp is the half width 
at zero intensity. Numbers in parentheses are lcr uncertainties in units 
of the last digits. c Line fitted to three components, with the center 
frequencies fixed to the calculated values of the three hyperfine com- 
ponents and the v cxp parameter fixed at 14.5 km s . 




Rest frequency (MHz) 

Fig. 2. HIFI spectra of IRC +10216 showing the first three rota- 
tional lines of H 35 C1 and H 37 C1 (black histograms; spectral res- 
olution is 1 .5 MHz for the 7 = 1 - and 3 MHz for the 7 = 2 - 1 
and 7 = 3-2 lines). Line profiles resulting from the radiative 
transfer model are shown in magenta. For the 7=1-0 line, the 
calculated intensity has been divided among the three hyperfine 
components according to their line strengths. 

collapse of the hyperfine structure and to the line widening, ex- 
pressed in frequency, with increasing 7. The 7=1-0 line of 
both isotopologues is composed of three hyperfine components 
but a good determination of the parameters of each component 
was not possible, so that we only give the total velocity inte- 
grated intensity in Table Q] The remaining lines of HC1 and that 
of HF have roughly parabolic shapes, which indicates optically 
thick emission unresolved by the beam of HIFI. The 7=1-0 
line of HC1 was observed in IRC +10216 by |Peng etal~| (2010) 
with the CSO. The smaller beam size of the CSO (13.5"), as 
compared with that of HIFI (33.9"), implies a line intensity en- 
hancement of a factor of (33. 9/13. 5) 2 ~6 for a point-like source. 
This is essentially the value observed for the 7=1-0 line of 
H 35 C1 after correcting for the beam efficiency of each telescope 
(0.5 for the CSO and 0.75 for HIFI at 625.9 GHz). Therefore, 
the emission size of the HC1 7=1-0 line must be substantially 
smaller than the 13.5" beam of the CSO. This is consistent with 
the source size expected for the HC1 7=1-0 line based on its 
critical density (>10 7 cm -3 ), which is only exceeded in regions 
inside 10 15 cm (0.5" at a distance of 130 pc). 



4. Radiative transfer modelling 

To learn more about the excitation conditions and abundance of 
HF and HC1 in IRC +10216, we performed radiative transfer 
calculations based on a multi-shell large velocity gradient (LVG) 
formalism. The physical model consists of a spherical envelope 
of gas and dust around an AGB star with a radius of 4 x 10 13 cm 
and an effective temperature of 2330 K. The adopted mass loss 
rate and distance to IRC +10216 are 2 x 10~ 5 Mo yr~' and 130 
pc, respectively. Further details about the physical model will be 
given in a forthcoming paper (Agundez et al. in preparation). 
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To model the emission of the HF and HC1 lines, we in- 
cluded the first 10 and 15 rotational states within the ground 
and first-excited vibrational states. The v = 1 — band lies 
at a wavelength of 2.5 /mi for HF and 3.5 //m for HC1. 
Level energies were computed from the laboratory spectro- 
scopic constants of HF (fOdashim a et al.l 1999; Ram et al. 1996) 
and HC1 dCazzoli & Puzzarinil 2004; IDe Natale et~aT1 1997; 
ILe Blanc et al.l 1994), neglecting the hyperfine structure for this 
latter species. Dipole moments used to compute the line inten- 
sities of HF are 1.826526 and 1.87368 Debye for pure rota- 
tional transitions within the v = and v = 1 states, respec- 
tively (|Muenter & Klemperer| 1970; IPiecuch et aTl 1998), and 



angular distance (") 

1 10 1 



W 10 



0.099735 Debye for the fundamental ro-vibrational transition 
dPine et al.l 1985). For HC1. we used 1.1 09 and 1.139 Debye for 
the v = and v = 1 states, respectively ( |De Leeuw & Dym anus 
1971: IKaise71 1970). while for the v = 1-0 vibrational band 
we ado pted 0.0729 61 Debye for H 35 C1 and 0.073049 Debye for 
H 37 C1 dPine et al.l 1985). State-to-state rate constants for rota- 
tional de-excitation of HF through inelastic collisions with para 
and ortho H2 were taken from lGuillon et al.l (201 1). whose cal- 
culations cover the first seven rotational levels of HF and ex- 
tend up to a temperature of 200 K. An ortho-to-para ratio of 3 
was adopted for H2. For HF-He collisions, we adopted the val- 
ues calculated bv lReese et al.l (2005). covering the first ten lev- 
els of HF up to 300 K. The use of collision rate constants with 
H2, which are larger than with He by up to one order of mag- 
nitude, is critical for the excitation analysis of HF. We did not 
extrapolate the rate constants in temperature, based on the small 
variation in the collision cross sections for temperatures above 
200-300 K (see lGuillon eTal~l 2008. 201 1). As rate constants for 
collisions of HC1 and H2, we adopted the values calculated by 
INeufeld & Greenl (1994) up to 300 K with He as collider, scaled 
up by a factor of 1.35 and without extrapolating in temperature. 
Collision rates of HC1 with H2 may differ from those with He, 
although it is difficult to quantify these differences. 

The computed line profiles are plotted in Fig. Q] for HF and 
in Fig.|2]for H 35 C1 and H 37 C1. The agreement between observed 
and calculated line profiles is very good, except for the J - 2- 
1 transition of H 35 C1 and of H 37 C1, which are predicted to be 
somewhat stronger than observed. The modest signal-to-noise 
ratio of these lines and the lack of collision rate constants for 
HC1 and H2 could explain the discrepancies. 

We derive an abundance relative to H2 in the inner regions 
of the envelope of 8 x 10~ 9 for HF and 8 x 10~ 8 for H 35 C1, 
with a H 35 C1/H 37 C1 abundance ratio of 3.3 + 0.3, consistent with 
that derived from previous observations of NaCl, KC1, and A1C1 
in IRC +10216 dCernicharo et al.l l987; ICernicharo et al.1 2000). 
The abundances of HF and HC1 (H 35 C1 + H 37 C1) are plotted in 
Fig. [3] as a function of radius, for a model in which the initial 
abundances is set by the observations, but its radial dependence 
is a model prediction. Both HF and HC1 are predicted to show 
a decrease in abundance in the outer layers, owing to photodis- 
sociation by interstellar ultraviolet photons. However, most of 
the contribution to the observed line intensities (~80 %) comes 
from circumstellar regions inside ~2 x 10 15 cm for HF, and 2- 
4 x 10 15 cm for HC1, such that only the abundances in regions 
within these radii are properly sampled by the observed lines 
(see light shaded region in Fig. [3). The excitation of the lev- 
els involved in the observed HF and HC1 lines is dominated by 
inelastic collisions with H2, with the lines being subthermally 
excited beyond ~2 x 10 14 cm where the gas density is below 10 9 
cm 4 . Infrared pumping via the v = 1 vibrational state plays a 
non-negligible but minor role, enhancing the line intensities by 5 
% (HF) and 15 % (HC1) with respect to the case where infrared 
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Fig. 3. Abundances relative to H2 of HF and HC1 (including 
H 35 C1 and H 37 C1) as a function of radius in IRC +10216. 
Abundances plotted as solid lines have a constant value in the in- 
ner envelope, derived from a fit to the observed lines (see Fig.Q] 
and Fig- [2j, and decrease in the outer envelope due to photodis- 
sociation, as computed with a photochemical model. The region 
of the envelope sampled by the observed lines of HF and HC1 
(see text) extends up to 3 x 10 15 cm (light shaded region). Dashed 
lines correspond to abundances computed by thermochemical 
equilibrium, which only holds up to ~ 3 R* (dark shaded region). 

pumping is neglected. Rotational lines within the v = 1 state 
are predicted to be too weak to be detectable with HIFI, with 
predicted antenna temperatures of only ~0.025 K for the v = 1 
7 = 3-2 transition of H 35 C1 and less than 0.005 K for all other 
such transitions. The abundances derived here are estimated to 
be uncertain by a factor of two for both HF and HC1. The er- 
ror of the observations being 10-20 %, 2-15 % due to the noise 
in the spectra (see Table [TJ and 10 % due to the calibration of 
HIFI, most of the uncertainty comes from the model. For HC1, 
the error could be significantly larger if the collision rate con- 
stants with H2 as collider are substantially different from those 
with He. An increase in the collision rate constants would lower 
the abundance derived for HC1. 



5. Discussion 

The presence of both HF and HC1 in the inner circumstellar 
regions of IRC +10216 is consistent with the thermochemical 
equilibrium (TE) calculations of the atmospheres of cool stars. 
In Fig. [3j we compare the abundances derived for HF and HC1 
with those computed under the assumption of TE, adopting solar 
elemental abundances ( Aspl und et al.| 2009). The assumption of 
TE is only valid from the stellar photosphere out to a radius of 
~3 R* (dark shaded region in Fig.[3}, beyond which the decrease 
in density and temperature causes the chemical reaction rates to 
drop rapidly and the molecular abundances to freeze out to the 
values of the TE region (Agundez & Cernicharo 2006). 

Hydrogen fluoride and A1F (whose abundance relative to 
H 2 is 7.5 x lQ- 9 ; |Agundez] 2009) are observed to be the main 
fluorine-bearing species in the inner envelope of IRC +10216. 
Thermochemical equilibrium calculations predict that HF is the 
major reservoir of fluorine in the TE region around both oxygen- 
and carbon-rich AGB stars, with the abundance of HF being es- 
sentially equal to the elemental abundance of fluorine. The abun- 
dance of HF in the atmospheres of AGB stars, derived from ob- 
servations of ro-vibrational lines in the near-infrared, has been 
used as a proxy of the abundance of fluorine (Jorisse n et al.l 
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1992). These authors found F abundance enhancements of up 
to a factor of 30 over the solar value, and interpreted this as ev- 
idence of 19 F stellar nucleosynthesis in AGB stars. The forma- 
tion of fluorine remains poorly understood with three different 
proposed sources: low-metallicity low-mass AGB stars, Wolf- 
Rayet stars, and neutrino spallation in core-collapse supernovae 
(see e.g. lLucatelloet al.1 201 D. IAbiaet a"Ll (2010) found that the 
analysis of the HF ro-vibrational lines made by Jorisser Tet al.l 
(1992) was affected by line blends with C2 and CN, such that the 
derived F abundances are only slightly enhanced (with a mean 
value of 0.2 dex) over the solar value, in closer agreement with 
theoretical predictions of low-mass AGB stellar models. In IRC 
+ 10216, we derive, through an alternative method (the obser- 
vation of the HF 7=1-0 rotational line), a fluorine abundance 
of 3.60 ± 0.30 (in the usual logarithmic scale with H having a 
reference value of 12), which is almost one order of magnitude 
lower than the solar value (4.56 ± 0.30; |Asplund et al.| 2009). 
This seems to indicate that in IRC +10216 there is an important 
degree of depletion of F onto dust grains, even larger than that 
found in interstellar clouds where the HF abundance relative to 
H 2 is (1.1 - 1.6) x 10~ 8 dSonnentrucker et al.l 2010; |Monje et al.| 
2011). Our ability to place constraints on the stellar nucleosyn- 
thesis of fluorine in IRC +10216 is hampered by the HF abun- 
dance we derive being a lower limit to the elemental abundance 
of fluorine. Nevertheless, the low HF abundance derived here 
does not point toward any F abundance enhancement over the 
solar value in IRC +10216, thus arguing in favour of the results 
of Ab ia et aLl (2010). Further observations of the 7=1-0 line of 
HF in other AGB stars will help us to understand the degree of 
fluorine depletion onto dust grains and any possible abundance 
enhancement due to stellar nucleosynthesis. 

In the case of hydrogen chloride, in the inner regions of IRC 
+ 10216's envelope the derived abundance, including both iso- 
topes, is 1.1 x 10~ 7 relative to H2, making it the major reser- 
voir of chlorine, followed by A1C1 whose abundance relative to 
H2 is 5 x 10~ 8 (Agundez 2009). Other chlorine-bearing species 
such as NaCl and KC1 have lower abundances, 1.3 x 10~ 9 and 
7 x 10~ 10 relative to H2, respectively (Agundez 2009). Together, 
these species have a total chlorine abundance of 8 x 10 8 rela- 
tive to the total number of H nuclei, i.e. one quarter of the solar 



abundance of CI (3.2 x 10~ 7 relative to H; |Asplund etaLl 2009) 



The abundance of chlorine should not be modified inside AGB 



stars as 35 C1 and il Cl are thought to be formed during explosive 



37 



oxygen burning in supernovae ( Woosley et al.| 1973). Therefore, 
the missing chlorine in the inner envelope of IRC +10216 may 
be in the form of neutral atomic CI, which together with HC1 and 
A1C1 is predicted to lock most of the chlorine according to TE 
calculations (see Fig. 2 of ICernicharo et aLl 2010b). or alterna- 
tively may have condensed onto dust grains. 



6. Concluding remarks 

The hydrogen halides HF and HC1 have been proven to be ubiq- 
uitous molecules in different types of interstellar environments. 
The observations presented in this Letter show that these two 
molecules are also important constituents of circumstellar en- 
velopes around evolved stars. They would be formed in the in- 
ner regions close to the star under thermochemical equilibrium, 
which predicts that HF and HC1 should be the major reservoirs 
of fluorine and chlorine, respectively, in the atmospheres of both 
O- and C-rich AGB stars. The abundances derived here for the 
carbon star envelope IRC +10216 are substantially lower than 
predicted by TE calculations, suggesting that F and CI are likely 
to be severely depleted onto dust grains, by ~90 % for HF and 



~75 % for HC1, although some chlorine may also be in the form 
of atomic CI. We also find that in IRC +10216 the abundance of 
F is probably not enhanced over the solar value by nucleosyn- 
thesis in the AGB star, although this conclusion is weakened 
somewhat by the HF abundance we derive being a lower limit 
to the elemental abundance of fluorine. Both HF and HC1 should 
be detectable with HIFI through low J pure rotational transitions 
in other AGB stars, especially in those with high mass loss rates, 
which have a large amount of emitting material and are difficult 
to observe in the near-infrared owing to the opacity of the dusty 
envelope. 
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